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Abstract 

Comparisons  of  computed  flowfield  predic¬ 
tions  to  experimental  data  for  a  cone-cylinder 
configuration  at  Mach  6.86  and  angles  of  attack 
of  0,  6.7,  and  14  degrees  are  presented.  The 
two  inviscid  computer  codes  used  were  EAGLE  and 
ZEUS,  which  represent  a  time-iterative,  shock 
capturing  approach,  and  a  space-marching,  shock 
fitting  approach,  respectively.  Specific 
comparisons  were  made  between  the  computed  and 
the  experimental  surface  pressure  distributions 
along  the  body  at  various  circumferential  loca¬ 
tions  and  around  the  body  at  axial  locations  on 
the  conical  nose  and  the  cylindrical  afterbody. 
Comparisons  of  the  aerodynamic  oharacteristlca  of 
normal  force  coefficient  and  the  location  of  the 
oenter-of-pressure  were  also  made.  The  inviscid 
calculations  compare  well  with  experimental  data 
where  viscous  effects  are  minimal. 


Introduction 

The  renewed  interest  in  sustained  hypersonic 
flight,  coupled  with  the  lack  of  ground-based 
test  facilities  in  this  flight  regime,  has 
focused  attention  on  computational  fluid  dynamics 
(CFD)  codes  for  hypersonic  aerodynamic  studies. 
Although  inviscid  codes  obviously  cannot  address 
aerothermodynamlc  problems,  the  design  engineer 
can  use  these  codes  for  quick  and  efficient 
studies  of  aerodynamic  characteristics.  The 
purpose  of  this  paper  is  to  compare  the  hyper¬ 
sonic  flowfield  predictions  of  the  EAGLE  code  to 
the  ZEUS  code  and  to  experimental  data  for  a 
cone-cylinder  configuration  at  Mach  6.86  and 
angles-of-attac k  of  0,  6.7,  and  14  degrees. 
EAGLE  was  developed  for  application  in  the  high 
subsonic  and  transonic  regimes.  Confidence  in 
its  ability  in  the  supersonic  and  hypersonic 
regimes  will  expand  its  usefulness  as  a  design 
tool.  The  ZEUS  code  was  selected  for  comparison 
because  its  validity  in  high  Mach  ranges  is  well 
documented . 

In  the  next  section  a  description  of  the 
configuration  and  the  experimental  data  will  be 
presented.  The  following  section  will  give  a 
brief  discussion  of  each  of  the  codes  used, 
followed  by  a  section  presenting  results  of  their 
usage.  The  last  section  will  discuss 
conclusions . 


Configuration  and  Experimental  Data 

The  cone-cylinder  configuration  is  shown  in 
Figure  1.  The  conical  nose  had  a  half-angle  of 
10  degrees,  and  the  cylindrical  afterbody  was  4 
diameters  long.  The  cone-cylinder  was  tested  in 
the  Langley  11-inch  hypersonic  tunnel,  which 
achieved  uniform  flow  for  model  testing  at  Mach 
6.86  in  a  5-inch  square  test  section.  The 
Reynold's  number  was  290,000  based  on  maximum 
diameter.  A  detailed  description  of  the  hyper¬ 
sonic  tunnel  can  be  found  in  References  1  and  2. 


Figure  1.  Cone-cylinder  configuration. 


The  model  was  designed  with  a  length  of  8 
inches  in  order  to  retain  the  model  completely 
within  the  uniform  flow  region  of  the  test 
section  during  the  high  angle-of-attaok  cases. 
The  model  was  equipped  with  ten  pressure  orifices 
installed  along  the  generatrix  of  the  body,  five 
on  the  conical  nose,  four  on  the  cylindrical 
afterbody,  and  one  on  the  base.  The  model  was 
machined  from  steel  and  had  a  polished  surface. 


The  pressure  distribution  over  the  model  was 
obtained  at  circumferential  positions  correspond¬ 
ing  to  0,  30,  60,  90,  120,  150,  and  180  degrees 
for  the  6.7  and  14  degree  angle-of-attaok  oascr, 
and  only  at  the  0  circumferential  position  for 
the  0  angle-of-attack  case.  Force  and  moment 
measurements  for  the  cone-cylinder  model  were 
obtained  at  Intervals  of  about  3  or  4  degrees 
throughout  the  angle-of-attack  range.  A  thorough 
description  of  the  experiment  can  be  found  in 
Reference  3. 
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Computer  Codes 


EAGLE 

Program  EAGLE  (Eglin  Arbitrary  Geometry 
ImpLiolt  Euler)  ia  a  raultiblock  grid  generation 
and  ateady-atate  flow  solver  system.  This  system 
combines  a  boundary  conforming  surface  generation 
scheme,  a  composite  block  structure  grid  genera¬ 
tion  scheme,  and  a  raultiblock,  implicit  Euler 
flow  solver  algorithm.  The  three  codes  are 
Intended  to  be  used  sequentially  from  the 
definition  of  the  configuration  under  study  to 
the  flow  solution  about  the  configuration.  Full 
documentation  on  each  code's  usage  can  be  found 
in  References  H-7.  Program  EAGLE  has  been  used 
as  a  design  aid  in  the  analysis  of  both  free- 
stream  and  interference  aerodynamic  character- 
istios  of  single  and  multiple  arbitrarily  shaped 
bodies,  finned  and  unfinned,  in  the  high 
subsonic,  and  transonic  Mach  ranges. 

The  computational  grid  built  for  the  cone- 
cylinder  configuration  consisted  of  2  similar 
three-dimensional  blocks.  Each  block  covered 
half  the  cone-cylinder  configuration  and  was 
dimensioned  at  75  x  25  x  37.  That  is,  75  lines 
along  the  body,  25  lines  between  the  body  surface 
and  the  outer  boundary,  and  37  lines  circumfer¬ 
entially  about  the  body.  A  cross-section  of  the 
grid  is  shown  in  Figure  2.  MACA  TR  1135  tables'^ 
were  used  to  calculate  the  shook  wave  location 
from  the  body,  and  doubling  that  distance  was 
used  to  determine  the  distance  to  the  outer 
boundary  of  the  grid.  To  obtain  as  smooth  a  grid 
as  possible,  the  elliptic  grid  generator  was 
employed  for  100  Iterations. 


Side  View 


End  View 


Figure  2.  EAGLE  co<le  grid  -  75  x  25  x  37 
each  half. 


An  Euler  solver  was  used  to  predict  the 
pressure  distributions  and  force  and  moment 
calculations  for  the  cone-cylinder  configuration. 
The  solver  uses  an  implicit,  two-pass  upwind 
scheme,  and  is  second  order  accurate  in  space. 
It  uses  the  flux-difference  splitting  scheme  of 
Roe  and  is  stable  for  a  wide  range  of  Courant 
numbers  for  steady-state  computations.  All  cases 
here  were  run  with  a  Courant  number  of  i».  The 
Euler  equations  are  written  in  strong  conserva¬ 
tion  law  form  in  order  to  capture  discontinuities 
in  the  flow,  such  as  shocks.  A  finite  volume 
formulation  is  used  in  order  to  achieve  total 
flexibility  with  regard  to  geometry.  Local  time- 
stepping  is  used  to  accelerate  convergence  for 
steady-state  problems.  Characteristic  variable 
boundary  conditions  are  used  by  the  solver  on  far 
field  boundaries  and  at  the  body  surface.  As 
noted  earlier,  the  outer  boundary  was  placed  far 
enough  from  the  body  to  Include  all  shocks  within 
the  solution  domain,  so  freestream  flow  condi¬ 
tions  were  imposed  at  the  outer  grid  surface,  and 
extrapolation  was  used  at  the  outflow  plane. 

The  Euler  code  was  run  for  1,500  Iterations, 
which  coincidentally  corresponded  to  a  reduction 
of  the  L2  norm  of  the  residual  by  three  orders  of 
magnitude  from  its  initial  value.  This  was 
deemed  sufficiently  converged  for  comparison 
purposes.  The  0  and  6.7  degree  angle-of-attack 
cases  were  run  without  any  problems.  For  the  lU 
degree  angle-of-attack  oases,  the  restart  file 
from  the  1,500  iteration  6.7  degree  case  was  used 
to  nudge  the  angle-of-attack  up  to  8  degrees, 
then  10,  12,  and  finally  14  degrees.  Each 
successive  Increase  was  run  500  iterations. 

ZEOS 


An  inviscid,  finite  volume,  spatial  marching 
flow  solver  was  also  used  to  predict  surface 
pressures  on  the  cone-cylinder  configuration. 
The  Zonal  Euler  Solver  (ZEUS)  code,  developed  by 
Wardlaw,  et  al,'^  was  selected  because  of  its 
efficiency  and  robustness.  The  ZEUS  code  is 
suitable  for  computations  where  the  flowfield  is 
entirely  supersonic  about  the  body;  and  any 
canards,  wings,  or  fins  are  thin  relative  to  the 
body  diameter.  The  code  applies  a  zonal  strategy 
which  divides  the  computational  domain  into  one 
or  more  zones,  each  of  which  can  be  mapped 
separately.  Zone  edges  are  defined  by  the  user 
to  coincide  with  the  body  and  bow  shock 


surfaces. 
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The  calculation  of  a  body's  flowfield  is 
initiated  near  the  nose  tip  and  marched  aft  in  a 
planar  fashion  along  the  body  axis  to  the  base. 
The  Initial  starting  flowfield  is  generated  by  an 
external  source.  For  the  geometry  presented  in 
this  paper,  a  conical  starting  solution  was 
employed . 

The  ZEUS  code  was  used  to  calculate  flow- 
field  pressures  and  aerodynamic  coefficients  for 
this  cone-cylinder  configuration  at  Mach  6.86  for 
angles-of-attack  of  0,  6.7,  and  14  degrees. 


2 
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Pressures  between  the  body  surface  and  the  bow 
shook  were  calculated  at  grid  cell  centers. 
Surface  pressures  were  integrated  to  provide  the 
aerodynamic  normal  force  coefficient,  and  to 
determine  the  location  of  the  oenter-of-pressure . 

Since  the  cone-cylinder  is  axisymmetrie  with 
no  fins,  only  a  single  computational  zone  was 
used.  Due  to  symmetry,  the  zone  was  established 
between  the  leeward  and  windward  axes  around  one 
half  of  the  body. 

A  36  X  18  grid  was  employed  across  this 
zone.  There  were  36  grid  cells  equally  spaced 
between  the  leeward  and  windward  body  meridians, 
and  13  clustered  grid  cells  between  the  body 
surface  and  the  bow  shock.  These  cells  were 
clustered  radially  so  that  grid  points  were 
concentrated  close  to  the  body. 


The  marching  step  size  along  the  body  axis 
was  automatically  varied  as  a  function  of  the 
geometry,  flow  conditions,  and  Courant  (CFL) 
factor.  A  default  value  of  0.9  was  used  for  the 
CFL  factor  in  each  calculation. 


Results 


The  results  to  be  presented  in  this  section 
will  be  comparisons  between  the  experimental  and 
the  calculated  axial  and  circumferential  surface 
pressure  distributions.  Additionally,  the  cal¬ 
culated  normal  force  coefficient  and  the  location 
of  the  center-of-pressure  will  be  compared  to  the 
experimental  values.  Finally,  pressure  contours 
from  EAGLE  predictions  will  be  shown. 

It  should  be  noted  that  there  was  a  dis¬ 
crepancy  Identified  in  the  experimental  report's 
conclusions  of  the  quality  of  the  pressure  data. 
A  slight  deviation  of  the  pressure  on  the  fore¬ 
part  of  the  conical  nose  for  the  experimental 
data  is  considered  to  be  the  result  of  documented 
small  surface  irregularities  which  were  incurred 
in  machining  operations,  and  whose  effect  is 
accentuated  at  the  tip  where  the  imperfections 
become  relatively  large  in  comparison  to  the 
local  radius.’”  This  deviation  accounts  for 
irregularities  in  the  experimental  data  on  the 
nose . 


Figure  3  presents  .1  comparison  of  the  exper¬ 
imental  and  calculated  surface  pressure  distri¬ 
bution  for  the  leewarJ  symmetry  plane  at  0  degree 
angle-of-attack.  The  calculations  compare  favor¬ 
ably  with  the  data.  Note  the  EAGLE  code  over- 
predlcts  the  pressure  drop  across  the  Prandtl- 
Meyer  expansion  at  the  nose-cy 1 inder  Juncture 
more  than  ZEUS  does.  This  is  evident  in  all  the 
EAGLE  and  ZEUS  prediction  comparisons  to  follow. 


Figure  3.  Experimental  and  numerical 
comparison,  Mach  6.86, 

0  degree  angle-of-attaek. 

Figures  4,  5,  and  6  show  the  surface  pres¬ 
sure  comparisons  for  the  6.7  degree  angle-of- 
attack  case  at  the  leeward,  90  degree,  and 
windward  me-’idians.  At  the  leeward  meridian 
(Figure  4)  the  agreement  with  both  EAGLE  and  ZEUS 
is  poor.  This  is  due  to  cross-flow  effects  in 
the  experiment  which  generate  a  viscous  layer  on 
the  leeward  side.  At  the  90  degree  (Figure  5) 
and  windward  (Figure  6)  meridians,  each  code's 
surface  pressure  predictions  show  good  agreement 
with  the  experimental  data.  Discrepancies  at  the 
nose  tip  in  the  experimental  data  may  be  attrib¬ 
uted  to  model  defects  noted  earlier. 


Figure  4.  Experimental  and  numerical 

comparison,  Mach  6.86,  oa 

6.7  degree  angle-of-attbck, 
leeward  meridian. 
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Figure  5.  Experimental  and  numerical 
comparison,  Mach  6.86, 

6.7  degree  angle-of-attack, 
90  degree  meridian. 
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Figure  6.  Experimental  and  numerical 
comparison,  Mach  6.86, 

6.7  degree  angle-of-attack, 
windward  meridian. 


Figures  7  and  3  show  circumferential  surface 
pressure  comparisons  for  the  6.7  degree  angle-of- 
attack  case  at  an  axial  location  on  the  nose  and 
on  the  cylindrical  afterbody.  The  axial  location 
on  the  nose  (Figure  7)  was  at  X/L  =  .231,  and 
both  codes  show  a  slight  deviation  from  experi¬ 
mental  data  at  the  leeward  side  { clrcuraferentlaX 
position  of  0  degrees).  Comparisons  improve  near 
the  windward  meridian.  The  axial  location  on  the 
cylindrical  afterbody  (Figure  8)  was  at 
X/L  =  .758,  and  both  EAGLE  and  ZEUS  predictions 
deviate  quite  a  bit  from  experimental  data  around 
the  leeward  meridian.  Once  again,  predictions  on 
the  body  toward  the  windward  meridian  improve  and 
compare  quite  well.  These  two  figures  graphi¬ 
cally  display  the  cross-flow  effects  in  the 
experiment  which  generated  a  viscous  layer  on  the 
leeward  side.  As  those  effects  became  less 


pronounced  circumferentially  towards  the  windward 
meridian  the  codes  predict  surface  pressure  much 
better . 
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Figure  7.  Experimental  and  numerical 
comparison,  Mach  6.86, 

6.7  degree  angle-of-attack, 
X/L  s  0.231. 
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Figure  8.  Experimental  and  numerical 
comparison,  Mach  6.86, 

6.7  degree  angle-of-attack, 
X/L  =  0.758. 


Figures  9,  10,  and  11  show  the  surface 
pressure  comparisons  for  the  19  degree  angle-of- 
attack  case  at  the  leeward,  90  degree,  and 
windward  meridians.  At  the  leeward  meridian 
(Figure  9),  EAGLE  and  ZEUS  predictions  under¬ 
estimate  the  surface  pressures  from  the 
experiment.  Experimental  pressures  are  higher 
because  the  viscous  layer  has  slowed  down  the 
flow  on  the  leeward  side,  and  the  invisoid  codes 
cannot  catch  this  phenomenon.  ZEUS  predictions 
after  the  nose-cylinder  Juncture  show  a  little 
ringing  which  is  more  than  likely  due  to  a  step 
size  that  is  too  large  for  this  region.  At  the 
90  degree  (Figure  10)  and  the  windward  (Figure 


A 
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11)  meridians,  the  codes  show  good  agreement  with 
experimental  data.  Both  EAGLE  and  ZEUS  experi¬ 
ence  oscillations  of  the  flow  solution  at  the 
nose  tip. 


Figure  9.  Experimental  and  numerical 
comparison,  Mach  6.86, 

19  degree  angle-of-attaok, 
leeward  meridian. 


Figure  10.  Experimental  and  numerical 
comparison,  Mach  6.86, 

19  degree  angle-of-attack, 
90  degree  meridian. 
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Figure  11.  Experimental  and  numerical 
comparison.  Mach  6.86, 

14  degree  angle-of-attack, 
windward  meridian. 


Figures  12  and  13  show  the  circumferential 
surface  pressure  comparisons  for  the  19  degree 
angle-of-attack  case  at  an  axial  location  on  the 
nose  and  on  the  cylindrical  afterbody.  The  axial 
location  on  the  nose  (Figure  12)  was  at 
X/L  s  .231.  and  both  codes  show  a  slight  devia¬ 
tion  on  the  leeward  side  Just  as  the  6.7  degree 
angle-of-attaok  case  did.  The  axial  location  on 
the  cylindrical  afterbody  (Figure  13)  was  at 
X/L  5  .758,  and  the  same  type  of  correlation  is 
shown.  As  stated  for  the  previous  case,  the 
discrepancies  on  the  leeward  side  are  attributed 
to  cross-flow  effects  in  the  experimental  data 
which  generate  a  viscous  layer  on  the  leeward 
aide. 


Figure  12.  Experimental  and  numerical 
comparison,  Mach  6.86, 

19  degree  angle-of-attaok, 
X/L  =  0.231 . 


Figure  13.  Experimental  and  numerical 
oompariaon,  Mach  6.86, 
lit  degree  angle-of-attaolc, 
X/L  e  0.758. 


Figures  H(  and  15  present  comparisons 
tjetueen  the  experimentally  and  numerically 
derived  normal  force  coefficient,  and  the  loca¬ 
tion  from  the  nose  tip  of  the  center-of-pressure 
as  a  function  of  angle-of-attack.  Both  EAGLE  and 
ZEUS  predictions  of  normal  force  coefficient  were 
larger  than  the  normal  force  coefficient  derived 
from  the  experiment.  This  is  not  surprising  con¬ 
sidering  the  surface  pressure  distribution  com¬ 
parisons  on  the  leeward  side  of  the  configuration 
where  viscous  effects  are  more  pronounced.  Both 
codes  consistently  underpred icted  the  pressure 
distribution  at  and  around  this  meridian,  thereby 
making  the  normal  force  coefficient  slightly 
greater  than  the  coefficient  determined  In  the 
experiment.  notice  also  that  the  codes  success¬ 
fully  predicted  the  nonlinearity  of  the  normal 
force  coefficient  curve.  Figure  1A  shows  that 
both  EAGLE  and  ZEUS  predict  the  center-of- 
pressure  location  farther  aft  than  the  experi¬ 
mentally  determined  location,  although  the 
predicted  values  are  within  Si  of  the  experi¬ 
mental  values.  The  accuracy  of  the  calculated 
location  of  the  center-of-pressure  Is  a  good 
indication  of  the  accuracy  of  the  calculated 
pressure  distribution.  Both  EAGLE  and  ZEUS  do  a 
very  good  Job  calculating  pressure  distribution. 


Figure  Normal  force  coefficient 

comparisons . 


Figure  15.  Center  of  pressure 

location  comparisons. 


Figures  16,  17,  and  18  are  static  pressure 
contour  plots  from  EAGLE  flowfield  predictions 
for  the  0,  6.7,  and  li)  degree  angle-of-attack 
cases.  The  attached  shock-wave  angles  from  these 
contours  were  compared  to  NACA  TR  1135  values  for 
conical  flow  and  showed  excellent  agreement. 
Prand tl-Meyer  expansions  at  the  nose-cylinder 
Junctures  are  clearly  visible.  A  denser  overall 
computational  grid  or  clustered  grid  lines  in 
those  regions  of  sharp  flowfield  gradients  would 
yield  even  better  resolution  of  these  flow 
features . 
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Conclusions 


Figure  l6.  Pressure  contour  plots 
from  EAGLE  code,  Mach 
6.36,  0  degree  angle- 
of-attack. 


Figure  17.  Pressure  contour  plots 
from  EAGLE  code,  Mach 
6.86,  6.7  degree  angle- 
of-attack . 


Figure  18.  Pressure  contour  plots 
from  EAGLE  code,  Mach 
6,86,  111  degree  angle- 
of-attack . 


Hypersonic  flowfield  predictions  have  been 
computed  for  a  cone-cylinder  configuration  using 
EAGLE  code  and  ZEUS  code  at  Mach  6.86  and  0,  6.7, 
and  11)  degrees  angle-of-attack .  EAGLE  code 
predictions  of  surface  pressure  distributions, 
normal  force  coefficient,  and  center-of-pressure 
location  were  very  close  to  ZEUS  predictions,  and 
both  codes  compared  favorably  to  the  experimental 
data  in  those  regions  where  viscous  effects  were 
minimal.  These  results  represent  a  first  step 
toward  validating  the  EAGLE  code  in  supersonic 
and  hypersonic  applications.  This  validation, 
coupled  with  EAGLE  code's  validation  in  the  high 
subsonic  and  transonic  flight  regimes®''^,  will 
greatly  expand  its  usefulness  as  a  design  tool. 
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